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Abstract 

Based on the analysis of both Borel and Finite-Energy QCD sum rules, the 
inclusion of finite mesonic widths leads to a dramatic effect on the predictions 
for the momentum dependence of the p — uj mixing matrix element. It is 
shown that the p — uj mixing matrix element traditionally discussed in the 
literature, has the same sign and similar magnitude in the space-like region 
as its on-shell value. This contrasts the zero-width result where the mixing 
matrix element is typically of opposite sign in the space-like region. 
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Recently, it has been argued that the p — oj mixing matrix element should change signif¬ 
icantly off-shell [p]-[TT[|. This off-shell behavior would drastically reduce the contribution of 
p — uj mixing to the charge symmetry breaking (CSB) nucleon-nucleon potential In 

a recent Letter []^, we have shown, in a model independent way, that the hnite widths of p 
and UJ mesons can give rise to a new source of momentum dependence for the p — uj mixing 
matrix element. The dependence arising due to the meson widths leads to a signihcant 
alteration of the result obtained in the zero-width approximation usually discussed in the 
literature. The origin of this momentum dependence lies in the difference between the p 
and UJ meson widths. In Ref. [|T^ we have quoted various values of the mixing parameter A 
extracted from QCD sum-rule analysis that did not include mesonic widths . This Letter 
will explore the effect of hnite mesonic widths on the determination of A. 

To this end we invoke both the Borel QCD sum rules (BSR) and Finite-Energy QCD sum 
rules (FESR). We will show that the large difference between the p and uj widths can have 
a dramatic effect on the value of A. While it is positive in the zero width limit, A becomes 
negative when the physical p and uj widths are included. Adopting this hnite-width sum- 
rule prediction for A, we hnd that the p — uj mixing matrix element in the space-like region, 
—1.0 GeV^ < < 0, has the same sign and a similar magnitude as its on-shell value. This 

contrasts the zero-width result where the mixing matrix element is typically of opposite sign 
in the space-like region. 

Consider the p and uj mixed propagator dehned by 0] 


Ti%{q)=il dW^-(0|TJ,^(a;)J-(0)|0) , 


( 1 ) 


where Jjl and are interpolating helds representing the p and uj mesons, respectively. It 
has been argued recently that there is no unique choice for these interpolating helds |]9|,pR|. 
Although this is certainly the case, an obvious and physically preferred choice is the standard 
vector currents. The reader is referred to Ref. [jT^ for more discussions on this point. 

It is common to hnd two factorizations of the Lorentz structure of Il.^{q^) rehecting 
current conservation 


from which the relation 


n£;(9t = - (9^ - , 

(2) 

- q^l qu) n 2 (g^), 

(3) 

U,{q^)=q^n,{q^), 

(4) 


follows trivially. There has been a great deal of interest surrounding the value of ni(g^) 
at = 0 PJ7|JT0[|. The direct electromagnetic contribution (p —7 —>• a;) to the mixed 
correlator is reasonably well understood and makes a numerically small contribution to 
the correlator. As such, this physics is usually subtracted from the correlator |Q, and the 
resulting correlator focuses on the mixing of p and uj mesons as well as higher resonances. 
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It is then trivial to demonstrate that ni(g^) = 0 at = 0. In the absence of intermediate 
state photons, the correlator describes the propagation of physical hadrons. As snch, 

does not have any poles at = 0 as there are no hadronic resonances below the 
threshold of 4m^. If ni(g^) is hnite at g^ = 0 then relation (^) indicates that n 2 (g^) has 
a pole at g^ = 0. This would be perfectly acceptable provided there is a zero in the factor 
preceding 112 (g^) in (H) at g^ = 0 to cancel out the pole. This zero is absent in (||). Hence in 
order to correctly describe the propagation of physical hadrons, n 2 (g^) cannot have a pole 
at g^ = 0. From it follows that ni(g^) = 0 at g^ = 0. 

The analytic structure of the mixed propagator allows us to write a dispersion relation 
of the form. 


Reni(g 2 ) 


P r°° ImHifs) , 
“ / , 7- 

71 JAml [S — q^) 


(5) 


which is valid to leading order in quark masses with the neglect of small electromagnetic 
interactions between quarks 0. It is clear that there are many excited states contributing 
to ni(g^) in addition to the ground state p and oj resonances. However, the p — oo mixing 
matrix element is traditionally dehned by saturating the dispersion integral @ with the 
ground state p and oj resonances only [|] 


Re ni(g2) 


0 {q^) 

(g2 — rn?p){q^ — m'^) 


( 6 ) 


Consequently, the result, ni(0) = 0, is spoiled, and hence 0(0) 7 ^ 0. 

Saturating the dispersion integral of Eq. by the p and oj resonances with the inclusion 


of hnite mesonic widths, one can express the p — oj mixing element dehned above as P3 


^r(g^) ^ [ <^(9^) ^ Ip{q^) Ipjrri^) 

Ovi'nP) \G{m^) m? — Iw{q‘^) 

where we have dehned 






ds - 

(s - g 2 ) 


rrip^^ rp^^(g) _ 

(s - 
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(7) 

( 8 ) 


GA) ^ 


(g2 - mlf + ml Tl{qA [(g^ - mlf + ml T^ (g^)] 


(g 2 -m 2 )(g 2 


'ml) 


rrip m^ Tp{q^) r^(g2) 


(9) 


with m 2 = (m 2 + m 2 )/ 2 , m^ = and . 


.pj!-, ...p,^ — ..up^^ r 2 ,^/4, and Tp^^ is the half width of the meson. 
The momentum dependence of the mesonic widths is indicated in equations (P) and (^. In 
particular, the widths vanish below the threshold of 4m2. For simplicity we shall approx¬ 
imate Fp,^(g2) as constants given by their physical values, Fp = 151.5 MeV and F^^ = 8.4 
MeV. Inclusion of momentum dependent widths has a negligible ehect on the mixing in the 
space-like region of interest . 
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Our task now is to determine the mixing parameter, A, using the QCD sum-rule ap¬ 
proach. QCD sum rules for p — uj mixing are obtained by studying the correlation func¬ 
tion of ( 0 ) with JP = — d'f^d)/2 and + dj^d)/6. The sum rules of 

= idfidu — dfiuff) are derived for 12n^‘^(g^) using a dispersion relation. The 

BSR can be written as 

where M is the Borel mass and Cq, Ci, C 2 and C 3 have been given in the literature. 


Q!e 


Cl = 


3 



~ 0 




C 2 



- mrf(l -Fy) 
2-^7 


{m)o , 


C 3 = 


224717 

81 


asK {qq)l 


287r 


tteK {qq) 


2 

0 ) 


( 11 ) 


where 7 = {dd)^/(^ u)q — 1 , {qq)o = {{uu)o +{dd)o)/2, and ae and are the electromagnetic 
and strong coupling constants, respectively. Here we have truncated the operator product 
expansion (OPE) at dimension six and kept only the terms considered in the literature. 

The studies of p —a; mixing via the QCD sum-rule approach have been made by Shifman, 
Vainshtein, and Zakharov 0 , and subsequently by Hatsuda et al. 0. In these works, the 
mesonic widths were ignored, and the spectral density Imn^‘^(s) has been characterized by 
poles at the masses of the p, cu, p' and uj' mesons plus a continuum model starting at an 
effective threshold Sq ||I4|,^: 


—Im n^‘^(s) 

Tl 


fp6{s- ml) - ml) + fp> S{s - mj,) 

-fuj'd{s-ml,) + CoO{s- So) . 


( 12 ) 


While the assumption of poles for lowest resonances may be reasonable for the determination 
of resonance parameters, here we are concerned with subtle cancelations of the nearby p and 
UJ strengths. To include the effect of hnite mesonic widths, we replace the 6 functions in 
Eq. (P!2D by normalized Briet-Wigner distributions for the meson widths of the form 

1 rUvTv 


TT (s — 771^)2 -I- ml Tl 
The BSR can then be expressed as 

fpip{M^) - /. /.(MQ + /,, J,,(MQ - f^, 1,{M^) 


(13) 


= c„Af- 




C2 C3 


where we have dehned 


R(MQ = - 


ds 




vr J4ml (s - ml)^ + ml T^ 


(14) 


(15) 
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The quantity of interest, A, is related to fp and through 


A = 


9p 9o. 
12 


A 

m] 


k 


mt 


m 


Q{m? 


( 16 ) 


where Q{rn?) is the on-shell value of the p — uj mixing with the electromagnetic contribution 
subtracted. Following the earlier works we take (?p/ 47 r = 2.4, = ?>gp, and Q{w?) = 

-O.OOSlGeV^. 

In Ref. [0, the authors used the BSR of (pR]), the hrst derivative sum rule as obtained by 
taking the derivative of ( pUD with respect to 1/M^, a FESR, and the on-shell constraint in 
their analysis. While the absence of strong continuum model contributions makes the p — uj 
mixing correlator a reasonable candidate for FESR analysis, we note that consideration of 
( 0 ) alone is sufficient to extract the parameters of interest. This latter approach sidesteps 


the serious problems associated with derivative sum rules |]15|,h6 


The on-shell constraint given in Refs. |1l4|,^ is valid for the zero-width case only. For the 
hnite-width case it becomes 


mt 


f = —t 

JuJ 9 




ipK: 

IJm^' 


fp + 


12 


m: 


0 (m^ 


5'pW 


(17) 


The BSR analysis is carried out following the techniques for determining the valid Borel 
region as introduced in Ref. |^. We limit the continuum model contributions to 50%, and 
maintain the contributions of the highest dimensional condensates in the operator product 
expansion (OPE) to less than 10% of the sum of terms in the OPE. The valid Borel regime 
throughout this work is 1.15 < M < 2.45 GeV. 

Uncertainties in the extracted parameters are determined via the Monte Garlo based 


uncertainty analysis introduced in Ref. |]^; uncertainties on QGD parameters such as con¬ 
densate values are utilized to provide uncertainties in the OPE. Here we use — 47 r^(gg)o = 
0.62 ± 0.05GeV^, k = 2 ± 1 and as/ 7 r(lGeV^) = 0.117 ± 0.014 |^. The average quark 


mass is chosen to satisfy the Gell-Mann-Oakes-Renner relation. The two new parameters 
encountered in (pR|) include mu/rrid = 0.50 ± 0.25 and 7 = 0.008 ± 0.002. The uncertainty 


in the quark mass ratio is taken from the particle data tables [|^ and is much larger than 
the uncertainty of ±0.04 assumed in Ref. [^. In turn, these OPE uncertainties provide 
distributions for the extracted parameters, from which uncertainty estimates are obtained. 
Since the resultant distributions are not Gaussian, we use the more robust measure of the 
median for the central value, and extract standard error estimates away from the median. 


1. Zero-Width Case 

To begin our analysis, we £x the continuum threshold, Sq; at the p(1700) resonance mass 
and optimize the £t in this region by adjusting the strength of the p' and uj' contributions. 
We utilize the constraint of (^), and extract A via (|T^). Table | summarizes the optimal £t 
parameters. These results are in reasonable agreement with Ref. |^. 

Of course, the 0 meson should be included in the parameterization of the spectral density 
II| . To keep the same number of parameters in the £t we replace the individual contributions 
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TABLE I. Borel sum rule predictions for the p — oo mixing parameters. 


Parameter 

Zero Widths 

Physical Widths 




Constrained 

No Constraint 

fp X 102 

U X 102 

1 qi +0.02 

1 95 

1 89^°'°^ 

1 q7+0.0i 
-0.02 

0 363+0.052 
'0-'Jua_Q g2g 

U.2485_Q QQQg 

0.28+°;27 

n OQ +0.20 
^‘^'^-0.11 

u X 10' 

— 

0.07 ±0.04 

-0.04 ±0.04 

-0.04 ±0.05 

/p/ X 102 

-0.11 ±0.20 

— 

— 

— 

L' X 102 

-0.15 ±0.23 

— 

— 

— 

SO (GeV2) 

2.89 

3.20 

6.00 

5.52 

A 

1 63+0-34 
-^•O'^-0.55 

9 11 “^0.37 
^•'*'^-0.61 

—0 64 

U.U‘1_0 75 

—0 62 


of p' and oo' by an effective strength fp^i at the average p' and oj' mass. Optimization of the 
three parameters fp, and sends the strength from the p', a;' region to zero. Hence, 
we discard the parameter in favor of optimizing the continuum threshold. Strength 
from this region is very small and has little effect on the optimal value of A. 

The continuum threshold is not precisely determined by the sum rule. We hnd = 
1.79 ± 1.50 GeV. This is due to a lack of excited state strength in the sum rules as governed 
by the leading terms of the OPE. To avoid the consideration of unphysical values for sq; 
we optimize Sq for the central values of the condensates, and then £x it for the Monte- 
Carlo sampling of the QCD parameters throughout the following. We note that fixing sq 
to the optimal value, as opposed to keeping this parameter as a search parameter in the 
uncertainty analysis, has a negligible effect on the central values and uncertainties of the 
other £t parameters. 

2. Physical-Width Case 

For the hnite physical-width case, we hrst consider a £t of /p, and Sq and determine 
fuj from the constraint of ([T7|) . The central result of this paper is that A changes sign due 
to the large hnite width of the p meson. The value of sq is larger than that for the zero 
width case, indicating a reduction of strength in the excited state regime, including the p', 
uj' region. 

Using the on-shell constraint effectively fixes one of the parameters we are trying to 
extract from the QCD sum rules. Since there could be a discrepancy between the on-shell 
value and the QCD parameters we have selected, it is imperative to relax this constraint. 
Since the p-meson width is much larger than the a;-meson width, there is now a considerable 
difference between their contributions to the sum rules. Thus, it is possible to resolve 
the individual contributions of p and uj mesons without the external constraint. The “No 
Constraint” column of Table | summarizes the results of this four parameter fit. While the 
absolute values of fp and are not as well determined, their correlations lead to a value 
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for A that is essentially unchanged from the constrained case and with only slightly larger 
uncertainty. The sum rule prediction for the on-shell value is 0(m^) = —0.0046GeV^ 
which agrees with the physical value of —0.0051 GeV^. 

In searching for correlations between the QGD parameters of the OPE and the p — uj 
mixing parameter A, only one interesting correlation arises. Fig. displays a scatter plot 
for A and the quark mass ratio mu/md where the constraint equation of ([ITD is used. If 
the constraint equation is relaxed for the hnite width case, the distribution remains largely 
unchanged. Finally, if we adopt the value for m^/m^ = 0.57 ± 0.04 estimated in and 
used in [0, we hnd A = —0.52 ± 0.32, using the constraint of (pG|). 

Since the mixed correlator has little excited state strength, we have also performed the 
FESR analysis. Results are shown in Table y. We see that the FESR results are in good 
agreement with the BSR results. To emphasize that the change in sign of A occurs due to 
the large difference of p and uj widths, we show the results for a hypothetical case of hnite 
but equal meson widths in column 5. In this case the results agree with the zero-width case. 

The implications of hnite mesonic widths on the dependence of 6'r(q'^)/6'r(m^) are 
illustrated in Fig. The ratio 0Y{q^)/Oi-irn?) is positive in the space-like region in 
contrast to the zero-width results. For the median value of A, the p — uj mixing matrix 
element in the space-like region has the same sign and similar magnitude as its on-shell 
value. A comprehensive analysis of the impact of the present result on the GSB potential 
and corresponding physical observables will be documented in a future publication . 


In the present paper, the p—uj mixing matrix element is dehned to include only the ground 
state p and uj. The contributions of excited states may also be included in a more complete 
dehnition. However, this will require the knowledge of the couplings of the interpolating 
helds to the excited states. The feasibility of reliably determining such couplings via QGD 
sum rules is currently under investigation. 


We would like to thank Ernie Henley, Jerry Miller, and Richard Woloshyn for helpful dis¬ 
cussions and comments. This work was supported by the Natural Sciences and Engineering 
Research Gouncil of Ganada and the US Department of Energy under grant DE-FG06- 
88ER40427. 
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FIG. 1. Scatter plot of A obtained from 1000 QCD parameter samples. The correlation with 
the quark mass ratio is displayed. The effect of including finite mesonic widths (A) is to shift down 
the zero-width results (o) and produce negative values for A. 


TABLE II. Finite-Energy sum rule predictions for p — u) mixing parameters. 


Parameter 


Zero Widths 

Finite Widths 

Physical Widths 

Equal Widths 

fp X 102 

1.87 

1.85 

0.304 

1.86 

U X 102 

1.91 

1.93 

0.248 

1.93 

u X 10" 

— 

-0.081 

0.021 

-0.073 

fp’ X 10^ 

-0.17 

— 

— 

— 

L' X 10^ 

-0.21 

— 

— 

— 

fp'u.' X 102 

— 

-0.0003 

-0.0037 

0.0016 

so (GeV2) 

3.0 

3.0 

3.0 

3.0 

A 

1.53 

2.13 

-0.67 

2.02 
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(GeV^) 

FIG. 2. Or{q^)/Or{m^) as a function of for A = —0.64^Q’y® obtained from BSR analysis 
with the on-shell constraint. The solid curve corresponds to the median value, and the two dotted 
curves illustrate the standard error. The two spikes in the time-like region correspond to the two 
poles of G{q^) [see Eq. @]. 


REFERENCES 

[1] T. Goldman, J. A. Henderson and A. W. Thomas, Few-Body Systems 12, (1992) 123; 
Mod. Phys. Lett. A7 (1992) 3037. 

[2] G. Krein, A. W. Thomas, and A. G. Williams, Phys. Lett. B317 (1993) 293. 

[3] J. Piekarewicz and A. G. Williams, Phys. Rev. G 47 (1993) R2461. 

[4] T. Hatsuda, E. M. Henley, T. Meissner and G. Krein, Phys. Rev. G 49 (1994) 452. 

[5] K. L. Mitchell, P. G. Tandy, G. D. Roberts and R. T. Gahill, Phys. Lett. B335 (1994) 
282. 

[6] H. B. O’Gonnell, B. G. Pearce, A. W. Thomas and A. G. Williams, Phys. Lett. B336 
(1994) 1; Phys. Lett. B354 (1995) 14. 

[7] R. Urech, Phys. Lett. B355 (1995) 308. 

[8] A. N. Mitra and K. G.Yang, Phys. Rev. G 51 (1995) 3404. 

[9] T. D. Gohen and G. A. Miller, Phys. Rev. G 52 (1995) 3428. 

[10] K. Maltman, Phys. Lett. B362 (1995) 11. 


9 

























[11] S. Gao, C.M. Shakin, and W. Sun, Phys. Rev. C 53 (1996) 1374. 

[12] M.J. Iqbal and J.A. Niskanen, Phys. Lett. B322 (1994) 7. 

[13] M.J. Iqbal, X. Jin, and D. B. Leinweber, Phys. Lett. B367 (1996) 45. 

[14] M. A. Shifman, A. 1.Vainshtein, and V. 1. Zakharov, Nucl. Phys. B147 (1979) 519. 

[15] X. Jin and D. B. Leinweber, Phys. Rev. C 52 (1995) 3344. 

[16] D. B. Leinweber, preprint UW-DOE/ER/40427-17-N95, |nucl-th/9510051| (1995). 

[17] D. B. Leinweber, Ann. Phys. 198 (1990) 203. 

[18] L. Montanet et al, Phys. Rev. D 50 (1994) 1173. 

[19] K. Maltman, Phys. Rev. D 53 (1996) 2563. 

[20] J. Gasser and H. Leutwyler, Nucl. Phys. B250 (1985) 465. 

[21] M. J. Iqbal, X. Jin, D. B. Leinweber, in preparation. 


10 





